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Convection is the mechanism by which energy is trans- 
ported through the outermost 30% of the Sun Q]. Solar 
turbulent convection is notoriously difficult to model across 
the entire convection zone where the density spans many or- 
ders of magnitude. In this issue of PNAS, Hanasoge et al. 
[2] employ recent helioseismic observations to derive stringent 
empirical constraints on the amplitude of large-scale convec- 
tive velocities in the solar interior. They report an upper limit 
that is far smaller than predicted by a popular hydrodynamic 
numerical simulation. 

Historically, great advances in our understanding of the 
solar interior have been due to helioseismology, the study of 
five-minute solar internal oscillations [3]. In the mid 1980s 
global-mode frequencies were used to measure the depth of 
the solar convective envelope at 0.71 solar radius, deeper than 
previous expectations based on underestimated opacities. An- 
other spectacular achievement was the inference of solar ro- 
tation as a function of radius and latitude. The bulk of the 
convective envelope rotates differentially, faster at the equator 
than at high latitudes. At the base of the convection zone is a 
zone of rotational shear, known as the tachocline, which now 
plays a central role in theories of the solar dynamo [4]. Despite 
valuable attempts, none of the above solar features were confi- 
dently predicted by models. Whenever helioseismology opens 
a new window into the solar interior, surprises are possible. 

The work of Hanasoge et al. [2] is perhaps the most no- 
table helioseismology result since the launch of the Helioseis- 
mic and Magnetic Imager (HMI) [5] onboard NASA's Solar 
Dynamics Observatory (SDO). HMI measures the motions 
on the solar surface caused by the random superposition of 
seismic waves excited by near-surface convection. Full-Sun 
Doppler velocity images are captured every 45 seconds by a 
16 million pixel camera. Hanasoge et al. used this unique 
combination of high resolution and full spatial coverage to 
carry out high-precision helioseismology of large-scale solar 
convection. 

To first order, convective flows do not affect the global- 
mode frequencies of solar oscillations. On the other hand, in 
time-distance helioseismology [6], wave travel times are lin- 
early sensitive to subsurface flows [7]. Time-distance helio- 
seismology uses spatio-temporal correlations of the random 
wave field to measure the travel times of solar acoustic waves 
between distant locations on the solar surface. The two- 
point correlation function contains fundamental information: 
it plays a role similar to a Green's function. One applica- 
tion of time-distance helioseismology has been the study of 
solar supergranulation|8]; This preferred scale of convection 
with typical length scale 35 Mm has remained a challenge for 
theory since its discovery fifty years ago [5]. 

Hanasoge et al. [2] measured east-west travel-time differ- 
ences using a deep-focusing, quadrant geometry, for a range 
of integration times, T, up to 96 hr. Travel times are dom- 
inated by stochastic noise due to the incoherent superposi- 
tion of solar seismic waves. According to Hanasoge et al., an 
upper limit on convective velocities at harmonic degree t is 
v$,£ < e ss Te/Ct. In this expression, e sg <C 1 is the signal-to- 
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Fig. 1. Comparison of kinetic energy spectra E^ of longitudinal solar velocities 
versus spherical harmonic degree £. The black curve above the gray area shows the 
observational upper limit from SDO/HMI helioseismology at radius 0.96-Rq and 
one day averaging 2 . The red curve is from surface velocities measured by tracking 
granules using SDO/HMI intensity images 18 . Notice the excess power at i ~ 120 
due to supergranulation. Results from snapshots of numerical simulations of con- 
vection sliced at 0.98/2© are given by the blue and light blue curves. The blue 
curve shows the spectrum from an ASH spherical-shell simulation 1121 and the light 
blue curve shows the spectrum from a stagger 16 near-surface radiative compressible 
simulation of size 96 X 96 X 20 Mm^. The horizontal black line and associated 
arrow shows a theoretical lower limit based on global dynamics arguments 1141 . as- 
suming mode equipartition over £ < 750. We defined E^ at radius r such that 
(v^)/2 = ^2i>0 Efi{£)/r, where is the longitudinal (prograde) component 
of the velocity field and the expectation value is approximated by a horizontal average 
(consistent with 13 in the limit of large £). 



noise ratio of the travel times (dominated by the contribution 
from supergranulation) and ti is the rms travel time at scale 
i only (dominated by stochastic noise). The quantity e ag may 
be estimated from the T-dependence of the variance of the 
travel times, since the variance of pure noise behaves like 1/T 
[10] , The measurement of e sg by Hanasoge et al. [2] is an 
advance compared to the earlier approach of Hanasoge et al. 
[11] where e sg = 1 was assumed. The calibration constant 
Ct used above to convert travel times into velocities (see [2J, 
Figure 5 of supplementary material) is the result of numerical 
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forward modelling and leads to an upper limit on the velocity 
at the target depth. 

The interpretation of travel-time measurements is a topic 
of current research, especially regarding the effect of time- 
dependent turbulent velocities. Hanasoge et al. assumed 
frozen convection to obtain the calibration constant Ci, which 
is only an approximation when the lifetime of convection is 
less than T. This may lead to an underestimation of the so- 
lar velocities. Also the vertical correlations of the convective 
velocities were ignored: Eddies of vertical sizes less than the 
first Fresnel zone cannot be detected. These two points may 
potentially affect the inferred upper limit on the convective 
velocities. Furthermore, experience with other experiments 
suggests that systematic effects, e.g. center to limb effects, 
could also play a role at the m/s level. 

It is enlightening to consider the helioseismology results 
in the context of existing models of convection. Hanasoge 
et al. 2] showed that in the range I < 60 their helioseismol- 
ogy upper limit for longitudinal convective velocities at radius 
O.92J?0 is orders of magnitude less than what is predicted by 
an Anelastic Spherical Harmonic (ASH) hydrodynamic sim- 
ulation _12;. This is a main concern of the authors. For the 
outside commentator, there is no clear way to reconcile this 
severe disagreement. 

Here we supplement the comparison with additional the- 
oretical, numerical, and observational constraints, which we 
have combined in the figure. We chose the kinetic energy 
density E$(£) to characterize the strength of longitudinal con- 
vective velocities [13] , which is a partial but useful description 
of convection. 

Miesch et al. Q3] recently obtained an interesting the- 
oretical lower limit of 30 ms" 1 for convective velocities at 
radius 0.95Rq for scales £ < 750. This calculation is based on 
the idea that the observed large-scale flows (differential rota- 
tion and meridional circulation) are maintained by convective 
Reynolds stresses. Partitioning the kinetic energy evenly over 
all modes £ < 750, we find E$ > 0.4 km 3 s~ 2 , which is well 
above the helioseismology upper limits at the lowest £ values 
(see figure). More work is needed to determine the £ depen- 
dence of this theoretical lower limit. 

The ASH simulation is truncated at radius O.98J?0, above 
which additional physics is needed, for example compressibil- 
ity and radiative transfer. Convection in the near-surface 
layers has been modelled with great success (judging from 
comparisons with surface observations |15| ) using fully com- 
pressible radiative simulations in local Cartesian simulation 
boxes. For example a recent stagger simulation [16] covers 
r > O.977?0, which overlaps in radius with the ASH simu- 
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lation. The kinetic energy spectra of the two simulations at 
r — 0.98i?Q roughly agree around I ~ 150. This suggests 
that ASH is capturing some of the general dynamics there, 
despite using simplified physics and missing the driving by 
strong cooling at the solar surface. 

SDO/HMI provides at least two other means than helio- 
seismology to observe surface flows. The first method is direct 
Doppler measurements of the line-of-sight component of veloc- 
ity [17] . The second method is based on tracking the motions 
of granules or other small features and provides both com- 
ponents of the horizontal velocities. In the Figure, we plot 
the granulation-tracking result from Roudier et al. [18] , who 
employed SDO/HMI intensity images. Notice that the ASH 
kinetic energies are above the granulation-tracking value for 
£ < 80. This is surprising since one would expect the convec- 
tive velocities to decrease in amplitude with increasing depth. 

A striking feature in the granulation-tracking curve is the 
excess kinetic energy at I ~ 120 due to supergranulation. 
Current simulations of convection are not ideal for modelling 
the supergranulation in detail; currently this scale is near 
the largest scale of the stagger code and the smallest scale 
of ASH. The helioseismology inferences from Hanasoge et al. 
stop short of the supergranular scale; the method of anal- 
ysis was not optimized for that purpose. Overall, a better 
observational coverage and theoretical understanding of the 
intermediate spatial scales would help connect the local and 
global scales of convection. The next generation of convection 
models are expected to cover the supergranulation range. 

Assuming that the helioseismology upper limit on convec- 
tive velocities from Hanasoge et al. [2] can be taken at face 
value, this will force a rethinking of the large-scale dynamics 
of the solar convective zone. One particular question is how to 
model very highly turbulent regimes, e.g. by including deep 
thermal plumes [19] . 

Any viable theory of convection ought to explain convec- 
tion in other stars. In this respect asteroseismology may play 
an important role. The observed amplitudes of oscillation in 
other Sun-like stars [20] contain information about the vigour 
of surface convection in these stars, which in turn will place 
constraints on stellar convection models. 
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